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ABSTRACT 

We used the Magellan adaptive optics system to image the 11 Myr substellar companion IRXS 
1609 B at the bluest wavelengths to date {z' and Yg). Comparison with synthetic spectra yields 
a higher temperature than previous studies of Teff = 2000 ± 100 K and significant dust extinction of 
Ay = 4.51q‘7 Mass estimates based on the DUSTY tracks gives 0.012-0.015 Mq, making the 

companion likely a low-mass brown dwarf surrounded by a dusty disk. Our study suggests that IRXS 
1609 B is one of the ^25% of Upper Scorpius low-mass members harboring disks, and it may have 
formed like a star and not a planet out at ^320 AU. 

Keywords: brown dwarfs — instrumentation: adaptive optics — planetary systems — planets and 
satellites: individual (IRXS J160929.1-210524 B) — stars: individual (IRXS J160929.1- 
210524) — stars: pre-main sequence 


1. INTRODUCTION 

Discovery of substellar companions at hundreds of AU 
from their host stars in direct imaging surveys has posed 
challeng es to classical format ion mechanisms like co re ac¬ 
cretio n (jPollack et al.l 1 19961 ) and disk instability (jBossI 
Il997[ ) . The ultra-wide separations make the core-growing 
timescale exceedingly long, and observers have not yet 
discovered such long-lived protoplanetary disks. Some 
alternatives have been proposed, including gravitat ional 
scattering to the current location (|Veras et al.ll2QQ9[) and 
in situ star-like formation. 

Detecting and characterizing circumsubstellar disks 
can place constraints on formation mechanisms. Many 
substellar (close to planetary mass) companions have 
been suggested to host their own disks (separate from 
the primary’s disk) based on 0/IR emission lines or ex¬ 
cess. T he 1.28 /im Pa 0 em ission line w as seen on GQ 
Lup B (jSeifahrt et al.l [200^. CT Cha B JSchmidt et alJ 
[20(M . GSG Q6214-QQ21QB (iBowler et al.ll2011[) . and FW 
Tau C (iBowler et al.ll2Q14l). Hu bble Spaee Teleseope ob¬ 
servations by IZhou et ahl (j2Q14f ) also showed that GSG 
06214-00210 B, GQ Lup B, and DH Tau B exhibit an 
optical excess at ^0.3-0.7 /im, imp lying a mass accre¬ 
tion rate of 10“^^-10“^ Mq yr“^. iKraus et al.l (j2014f ) 
also found that IRXS 1609 B, GQ Lup B, GSG 06214- 
00210 B, DH Tau B, and ROXs 12 B have redder K' — V 
colors than young field dwarfs. An unresolved 24 /im ex¬ 
cess was also detected on G SG 06214-00210 B and IRXS 
1609 B (jBailev et al.ll201^ . In particular, iKraus et al.l 

* As argued bv IPecaut et al.l (|20I2I ). GSC 062I3-03I58 B or 
[PZ99] J160930.3-210459 B might be better names. 

^ This paper includes data gathered with the 6.5 m Magellan 
Clay Telescope at Las Campanas Observatory, Chile. 

® NASA Sagan Fellow. 

^ Carnegie Postdoctoral Fellow. 


(j2015f ) presented an ALMA Gycle 1 1.3 mm detection of 
dust continuum emission associate d with FW Tau C and 
derived a dust mass of 1-2 Mq. iGaceres et al.l (|2015f ) 
further detected GO (2-1) and showed that gas is also 
present in this disk. Although FW Tau G could be in 
the b rown dwarf regime (M ^ 10 ± 4 Mjup; IKraus et al.l 
EoH), these observations may represent the first direct 
detection of a disk around a planetary mass object. 

Recently, usi ng the 6.5m Mage ll an adaptive optics sys¬ 
tem (MagAO; iClose et~al] 120131: iMales et al.l I2014D we 
have also detected an r' (0.63 ptm) excess possibly due 
to Hg and a dust extin ction of Ay = 3-4 mag for GT 
Gha B (jWu et al.l 120151 ). All of these observations sug¬ 
gest that circumsubstellar disks could be common around 
these wide young planetary mass objects. The existence 
of disks also favors star-like fragmentation, but argues 
against a scattering origin. This is because disks may 
be perturbed, if not de stroyed, by each encounter with 
anoth er massive body (jBowler et al.l [20TTI: iBailev et al.l 

iMi. 

Here we present Mag AO z^ and Yg imaging of IRXS 
1609 B, a substellar companion discovered at 320 AU 
(projected separation) from IRXS J 160929.1-210524 
in the Upper Scorpius association by iLafreniere et al.l 
(1200^ ). The companion is widely reported as the 
first directly imaged exoplanet orbiting a Sun-like star. 
Its mass, temperature, and spectral type were deter- 
mined to be 0.008-0 .0 11 M ^ d, ^1800 K, and ^L 4 
(jLafreniere et al.l I2008L 120101: iLachapelle et al.l 120151) . 
Recently, iPecaut et al.l (j2012f ) revised the age for Upper 
Sco to be 11 ± 2 Myr, and, based on that, they derived 
a higher mass 14^3 Mj^p. In this paper, we show that 
IRXS 1609 B may have some dust extinction, a higher 
temperature, and an earlier spectral type. Therefore, it 
is more likely a low-mass brown dwarf slightly obscured 
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by a dusty circumsecondary disk. 

2. OBSERVATIONS AND REDUCTION 

We used Mag AO, a new AO system on the 6.5 m Clay 
Telescope, to image IRXS 1609 A and B at z' (A = 
0.91 /im; AA = 0.12 fim) and Yg (A = 0.98 /im; 
AA = 0.09 jam) on 2013 April 6 (UT) during the sec¬ 
ond commissioning run. Seeing varied between 0'.'6 and 
0'.'8. We used the primary {R ^ 12.4 mag) as the guide 
star, and locked the AO system at 100 Karhunen-Loeve 
modes and 400 Hz. This is a very faint target given 
that 50% of the light goes to the VisAO science cam¬ 
era and 50% to the pyramid wavefront sensor. That is 
why only 100 of 378 possible modes were corrected at 
400 Hz compared to the usual 1000 Hz loop speed. The 
resulting FWHMs were 67 and 72 mas for z' and Tg, 
respectively. For z', we obtained 20 s x 193 (3860 s) 
and 2.27 s x 28 (63.6 s) for saturated and unsaturated 
data, respectively. For Fg, we obtained 20 s x 126 (2520 
s) unsaturated frames. We only detected IRXS 1609 B 
at z' but not Fg, possibly due to a lower quantum ef¬ 
ficiency which leads to a lower total throughput (0.013 
versus 0.077). To calibrate our photometry, we retrieved 
K7V, MOV, and MIV spectral templates from the Pickles 
Atlas (|Picklesl 1199811. reddened them by Ay = 0.1 mag 
(extincti on of ATlBowler et al.ll20f3l . and integrated the 
DENIS (|EDchtein et al.|[l997 1 / and our z' and Yg filter 
curves over these templates as well as the Vega spectrum 
to derive I — z' and I — Yg colors. Then we applied these 
colors to the existing DENIS I measurement on the pri¬ 
mary to obtain z' and F^ photometry. We adopted 0.06 
and 0.10 mag as the uncertainties for z' and F^ based 
on comparison to E7V and MIV templates. We also ob¬ 
served the optical standard star LTT 3864 in the same 
observing run, but difficulties in dealing with Strehl ratio 
variation on different nights and the inherent noisiness of 
large aperture CCD photometry made calibrations much 
less precise than using the DENIS photometry. However, 
we found consistent results to within 10% for the pho¬ 
tometry for IRXS 1609 A. This demonstrated that the 
primary is not a very active variable and that the DENIS 
/ band photometry is valid for the night of 2013 April 6. 

In the following analysis on the z' data, we selected 
frames with wavefront errors less than 175 nm rms, cor¬ 
responding to the best two-thirds of data (129 fra mes). 
Data reduction were as detailed in IWu et aH (|2015[) . We 
constructed a master point spread function (PSE) from 
the unsaturated images and performed PSE-fitting pho¬ 
tometry. In addition, aperture photometry was done 
using an aperture of 1 EWHM in radius. Our final z' 
flux is the average of both approaches, and its uncer¬ 
tainty included a ^0.1 mag offset between them. To 
estimate any possible flux loss in halo subtraction, we 
subtracted scaled-down PSEs at the position of the com¬ 
panion (negative PSE injection) without removing the 
halo and found a consistent Az' with our aperture and 
PSE-fitting photometry. Therefore, we concluded that 
there was no significant flux loss when removing the ra¬ 
dially symmetric PSE halo profile. The astrometric erro r 
budget also included image distortion (|Wu et al.ll2Q15[) . 
Table [T] summarizes the system properties. 

3. RESULTS 

3.1. Properties of the Companion 


Table 1 

Properties of IRXS 1609 System 


Property 

Primary 


Companion 

Distance (pc)^’^ 


145 ± 14 


Separation ("p 


2.21 ± 0.01 


PA (°)^ 


27.1 ±0.3 


Age (Myr)*^ 
Spectral type 

MOT I® 

11 ±2 

L2±l^ 

Teff (K) 

^udu_2oo 


2000 ± 100 ^ 

Ay (mag) 

n 1 +0.3e 


4 pc + O.Sc 

log(Lbol/Lo) 

-O.STitO.lSf 


-3.36 ± 0.09^ 

Mass (Mq) 

0.85l°;20t 


0.012-0.015^ 

IS 

10.99 ±0.03 



z'^ 

10.60 ±0.06 


21.24 ±0.15 

w 

10.43 ±0.10 


>19.46 (3o-) 

Jh,i 

9.764 ±0.027 


17.85 ±0.12 


9.109 ±0.023 


16.86 ±0.07 


8.891 ±0.021 


16.15 ±0.05 

3.1 /rm) 

8.80 ±0.05 


15.65 ±0.21 

3.3 /j,m^ 

8.78 ±0.05 


15.20 ±0.16 

L'k 

24 iim (mJy)j 

8.73 ±0.05 

3.06 ±0.04 

14.8 ±0.3 


Note. — _ _ _ 

^ Ide Zeeuw et al.l (|I9990 . ^ llreland et al.l (120111) . ^ This 

work. ^ IPecaut et al.l (| 20 I 2 l b ® IBowler et al.l (|20I4l b ^ 

ILafreniere et al.l (|2008l) . ^ DENIS (jEDchtein et al.lI997l) . ^ 2MASS 
(jSkrutskie et al.l 120061) . ^ ILachaoelle et al.l (|20I5I ) . ^ IBailev et al.l 
(|20I3I) . ^ ILafreniere et al.l (| 20 I 0 I ) . 



Figure 1. Left: IRXS 1609 in MagAO z' filter. Right: after 
subtracting the radially symmetric profile of the primary star. 


3.1.1. Temperature and Extinction 

Eigure [T] shows MagAO z' images of the system. 
We find a contrast of /S.z' = 10.64 ± 0.14 mag be¬ 
tween components. In Eigure O we compare our 
z' and Yg measurements, togeth er with the published 
JHKg spectra and photome try (|Lafreniere et al.l [2008L 
120101: iLachapclle et al.l 1201511. t o the 1 800 K Phoenix- 
ACES (jBarman et al ] \WrtR 1201 m . AMES-Dusty 
( Allard et, al.ll200lh . and four L7 dwarf^ in lAllers fc Liul 

(loil): 2MASS J05184616-2756457 (LI7), 2MASS 
J05361998-1920396 (L27), 2MASS J22081363+2921215 
(L37), and 2MASS J05012406-0010452 (L47). While the 
models and these Ly spectra fit IRXS 1609 B reason¬ 
ably well in the near-infrared, they all seem to be too 
bright at z' by a factor of ^2-4, suggesting that some 
dust might be present to redden the companion. There- 

® ICrnz etTin (120091) proposed a classification scheme of 7 and 
/3 for ~I0 Myr (low surface gravity) and ~I00 Myr (intermediate 
surface gravity) dwarfs, respectively. 
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Figure 2. Left: we compare the photometry and spectra of IRXS 1609 B (jLafreniere et al.l l2008U 2010l ') to the 1800 K AMES-Dusty and 
Phoenix-ACES models (Ay = 0 mag) and four low-surface gravity Ly dwarfs in lAllers' fc Liul (|2013IL We normalize them at Kg. Squares 
represent synthetic fluxes at z'YgJHKg^ and the Yg (0.98 /j,m) arrow extends from 3cr down to Icr. The companion’s z' flux is lower than 
that of Ly dwarfs and the models as well, suggesting that some dust extinction might be needed. We also note that these Ly objects may 
themselves also be reddened by dust; for example, the L2y appears redder than the L3 and L4 counterparts. 

Right: SED fitting with the 2000 K models. The blue curves represent what the object’s spectra would be if there was no extinction. 
Circles denote model fluxes for Alters. The Yg arrow extends from 3 f 7 down to Icr. We need a high extinction Ay ~ 4.5 mag (red curves) 
to match our observed Yg upper limit and z' flux. Note that a weak (~20%) ~3 /rm excess becomes apparent when the extincted SEDs 
are compared to the observations (3.05, 3.1, and 3.3 /xm data points are slightly above the red curves). 



J-H 



Figure 3. z — J vs. J — H of IRXS 1609 B a nd held L, 
T dwarfs (|Golimowski et ^120041 : IKnaop et al.ll 200 ^ : I Chiu et al.l 
1200611 . Once de-reddened by Ay ~ 4.5 mag, the color of IRXS 
1609 B is similar to that of L dwarfs. 

fore, in this study we explore other possibilities. We 
carry out a spectral energy distribution (SED) fitting us¬ 
ing three atmospheri c models: Phoenix -ACES, AMES- 
Dusty, and BT-Settl (|Allard et al.l 120111).^ We adopt log 
q = 4 .0 from previous analy ses (|Lafreniere et al.l [2QQ8L 
I2Q1QI: iLachapelle et al.l 12015!) . Then we ap ply the ex¬ 
tinctio n law in I Weingartner fc Draind (|2QQlf) an diD rainel 
( 2003 ) to redden these synthetic spectra. As in IWu et al. 
(|2015f) . to evaluate the goodness of fit, we match these 
reddened models with the observed Kg flux because the 
K band is most likely dominated by photospheric emis¬ 
sion and least affected by dust emission and extinc¬ 
tion. In general, we find that the reddened Phoenix- 


Figure 4. H-R diagram with DUSTY evolutionary tracks. Iso¬ 
mass curves are shown in gray. The salmon colored polygon is 
the only region consistent with the observed Tgff = 1900-2100 K, 
log(Lboi/L 0 ) = —3.36 ± 0.09, and the age of Upper Sco (11 ± 
2 Myr). This yields a mass estimate of 0.012-0.015 Mq. 

ACES synthetic spectra give the best fit, followed by the 
BT-Settl and AMES-Dusty models. Chi-square analysis 
gives Teff = 2000 ± lOOK, slightly higher than previous 
estimates. Extinction is more model dependent due to 
different treatments of opacity, varying from Ay = 2.7 
to 5.3 mag at this temperature range, so we compute a 
weighted average and adopt Ay = 4.5^0 ‘7 mag. 

Figure [2] also shows our SED fitting with models added 
Ay = 0 and 4.5 mag. We note that these fluxes are not 
absolute in order to avoid the ^10% distance uncertainty. 
We also normalize the companion’s JHKg spectra (gray 
curves) to their apparent fluxes. The red curves are nor- 
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malized at the apparent Ks flux, while the blue curves 
represent the object’s true SED when extinction is not 
present. Overall these reddened models fit the photome¬ 
try and spectra reasonably well. The 3 cr Yg upper limit 
also suggests that IRXS 1609 B is likely obscured, oth¬ 
erwise we would have detected it. 

We also notice that th e system has an u nresolved 24 
/im excess of 0.91 mJy (jBailev et al.ll 201 ^ . Since the 
primary star has v ery little extinction. Ay ^ 0.1 mag 
(jBowler et al.ll20l3 ) , it is possible that most of this warm 
dust excess originates from the compani on. _ As men¬ 

tioned in the Introduction, iKrans et al.l (j2014l ) found 
that IRXS 1609 B and a few other objects exhibit K' — L' 
excess, which indicates that disks may be common to 
young substellar companions. In Figure [ 2 ] there is evi¬ 
dence (<2 cr) of a weak ^20% 3-4 /im excess compared 
with our best-fit SEDs. In addition, in Figure [3] we 
find that I RXS 1609 B has 2 : — J redder than field L 
dwarfs (e.g . . iGolimowski et al.ll2003: iKnapp et al.l 12003 
iChiu et ahl |2006[) . Finally, the SED appears to need 
pure ISM-like reddening {Ay ^ 4.5 mag), in contrast 
to the gray extinction usually invoked in atmospheric 
planetary models wi th thick cloud s (e.g., HR 8799 b, 
Barman et al. 12011^ 2 M 1207 b, iSkemer et al.l I 201 TI: 
Barman et al.l^llbl ). All of these lines of evidence seem 
to suggest that IRXS 1609 B hosts its own inclined dust 
disk. Really, this is not a surprising result given that 
^25% of low-mass objects in Upper Sco have dusty disks 
(jLnhman fc Mama ield 120121 ) and IRXS 1609 B was the 
faintest, reddest o bject found by the disc overers’ AO sur¬ 
vey of Upper Sco (jLafreniere et al.ll2014l ). so the odds are 
good that the faintest, reddest object might also have 
dust extinction as well — making IRXS 1609 B appear 
redder and fainter than its intrinsic true color and lumi¬ 
nosity. 

3.1.2. Spectral Type, Luminosity, Mass, and Radius 

We calculate four gravity-insensitive i ndices H 2 O, 
H 2 OD, H 2 O-I, and H 2 O -2 as defined in lAllers fc Liul 
(j2013[ ). The average of these indices shows ^L 2 for Ay 
between 3 and 6 mag, earlier than previous estimates 
of L4 ± 1 (jLafreniere et al.ll2008L I 2 OIOI: iLachapelle et al.l 
l2015f ). Our result is also in a g reeme nt with the recent 
analysis by iManjavacas et all (j2014f ). who fit an L 27 
spectra to that of IRXS 1609 B. Therefore, we adopt 
L2±l. 

With Ay = 4 . 5 lo ‘7 D = 145 ± 14 pc, we cal¬ 

culate \o^(Lh a]/Lc^) = —3.36TO.09 using the bolometric 
correction in fSchmidt et al ] (j2014f ). To validate, we in¬ 
tegrate the de-reddened spectra in Figure [ 2 ] and obtain 
^ —3.32, suggesting that the bolometric correction de¬ 
rived from field dwarfs works reasonably well for young 
objects. Compyed t o the DUSTY evolutionary tracks 
(jChabrier et al.l[ 2000 l: see Figure E]), the companion has 
a mass betw e en 0.0 12 and 0.015 Mq, consistent with 
iPecant et al.l (j 2012 f ) but higher than 0.008-0.011 Mq 
found in other studies. Therefore, we suggest that IRXS 
1609 B likely lies above the fiducial brown dwarf/planet 
boundary. 

Finally, we obtain ^1.7 Jupiter radii using the new 
luminosity and temperature, consistent with the DUSTY 
tracks. 

3.2. Implications 


If IRXS 1609 B harbors an inclined disk, this will imply 
that circumsubstellar disks could survive after 10 Myr. 
This is not entirely unexpected because the recent large 
infrared survey in the Upp er Sco revealed longer dis k 
lifetimes for low-mass stars (jLnhman fc Mamajekll20f^ . 
The survival of disks also supports in situ fragmentation 
for companions on wide orbits and disfavors the planet- 
planet scattering scenario. Since no accretion-indicating 
lines were detected in the NIR spectrum, ongoing accre¬ 
tion is either slow or non-existent. The putative disk 
may be largely gas depleted, precluding accretion, while 
still retaining sufficient dust mass at larger radii to pro¬ 
duce the observed extinction. Future ALMA observa¬ 
tions could definitively test the existence of this disk. 
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